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On the Dilatational Stability of the Earth . 

By Lord Rayleigh, O.M., Pres. R.S. 

(Received March 14,—Read March 29,—Revised April 16, 1906.) 

The theory of elastic solids usually proceeds upon the assumption that the 
body is initially in a state of ease, free from stress and strain. Displacements 
from this condition, due to given forces, or vibrations about it, are then 
investigated, and they are subject to the limitation that Hooke's law shall 
be applicable throughout and that the strain shall everywhere be small. 
When we come to the case of the earth, supposed to be displaced from 
a state of ease by the mutual gravitation of its parts, these limits are trans¬ 
gressed; and several writers* who have adopted this point of view have 
indicated the obstacles which inevitably present themselves. In his 
interesting paperf Professor Jeans, in order to attain mathematical definite¬ 
ness, goes the length of introducing forces to counteract the self-gravitation: 
“ That is to say, we must artificially annul gravitation in the equilibrium 
configuration, so that this equilibrium configuration may be completely 
unstressed, and each element of matter be in its normal state." How wide 
a departure from actuality is here implied will be understood if we reflect that 
under such forces the interior of the earth would probably be as mobile 
as water. 

It appears to me that a satisfactory treatment of these problems must 
start from the condition of the earth as actually stressed by its self¬ 
gravitation, and that the difficulties to be faced in following such a course 
may not be so great as has been supposed. The stress, which is so enormous 
as to transcend all ordinary experience, is of the nature of a purely hydro¬ 
static pressure, and as to this surely there can be no serious difficulty. After 
great compression the response to further compressing stress is admittedly 
less than at first, but there is no reason to doubt that the reaction is purely 
elastic and that the material preserves its integrity. At this point it may be 
well to remark, in passing, upon the confusion often met with in geological 
and engineering writings arising from the failure to distinguish between 
a one-dimensional and a three-dimensional, or hydrostatic, pressure. When 
rock or cast iron is said to be crushed by such and such a pressure, it is the 

* See, for example, Love, 4 Theory of Elasticity/ § 127 ; Chree, 4 Phil. Mag., vol. 32, 
p. 233, 1891 ; Jeans, 4 Phil. Trans., 5 A, vol. 201, p. 157, 1903. 

t Loc. cit., p; 161. 
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former kind of pressure which is, or ought to be, meant. There is no evidence 
of crushing under purely hydrostatic pressure, however great. 

'Not only is the integrity of a body unimpaired by hydrostatic pressure, hut 
there is reason to think that the superaddition of such a pressure may 
preserve a body from rupture under stresses that would otherwise he fatal. 
FitzGerald raises this question in a review* of Hertz's c Miscellaneous Papers.' 
He writes : “In considering the cracking of a material like glass, Hertz seems 
to think its cracking will depend only on the tension; that it will crack 
where the tension exceeds a certain limit. He does not seem to consider 
whether it might not crack by shearing with hardly any tension. It is 
doubtful whether a material in which there were sufficient general compres¬ 
sion to prevent any tension at all, would crack. Eocks seem capable of being 
bent about and distorted to almost any extent without cracking, and this 
might very well be expected if they were at a sufficient depth under other 
rocks to prevent their parts being under tension. It is an interesting 
question whether a piece of glass could be bent without breaking if it were 
strained at the bottom of a sufficiently deep ocean. On the other hand, 
there seems very little doubt that the parts of a body might slide past one 
another under the action of a shear, and would certainly crack unless there 
were a sufficiently great compressional stress to prevent the crack ; and that 
consequently a body might crack, even though the tensions were not by 
themselves sufficiently great to cause separation, and might crack where the 
shear was greatest, and not where the tensions were greatest." 

When we reflect that pieces of lead may be made to unite under pressure 
when the surfaces are clean; and upon what is implied when insufficiently 
lubricated journals, or slabs of glass under polish, seize , we may well doubt 
whether it is possible to disintegrate a material at all when subjected to 
enormous hydrostatic pressure. In the words of Dr. Chreef: “ The con¬ 
ditions under which the deep-seated materials of the earth exist are 
fundamentally different from those we are familiar with at the surface. 
The enormous pressure, and the presumably high temperature, very likely 
combine to produce a state to which the terms solid, viscous, liquid, as we 
understand them, are alike inapplicable." 

A study of the mechanical operations of coining and of stamping (in 
recent years, I believe, much developed) would probably throw light upon 
this question. We know that rod or tube may be “ squirted " from hot (but 
solid) lead. Is the obstacle to a similar treatment of harder material purely 
practical? In the laboratory I have experimented upon jellies of various 

* i Nature,’ November 5, 1896 ; ‘ Scientific Writings/ p. 433. 

t ‘Phil. Mag./ vol. 43, p. 173, 1897. 
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degrees of stiffness, on the principle of suiting the material to the appliances 
rather than the appliances to the material. In the simplest arrangement 
a leaden bullet is imbedded in jelly contained in a strong glass tube which 
the bullet somewhat nearly fits. Although the tube stand vertical for 
several days, there is no appreciable descent. But if by numerous longi¬ 
tudinal impacts against a suitable pad the inertia of the bullet be brought 
into play, movements through several inches may be obtained. Here, 
although the deformations are very violent, there is no rupture visible, either 
before or behind the bullet. 

When an elastic body is slightly displaced from the condition of ease, the 
potential energy (V) is expressed by terms involving the squares and products 
of the displacements. If, however, we suppose given finite forces to be 
constantly imposed, so that the initial condition is one of strain, the case 
is somewhat, though not essentially, altered. It may be convenient to make 
a statement, once for all, in terms of generalised co-ordinates. If under the 
action of the forces <£>,©,... the co-ordinates assume the values <j>, 0,... we 
have in terms of the potential energy of strain Y, 


<$> 


etc. 


a) 


dV = dV 

d<j>’ d0‘ 

If the forces permanently imposed be distinguished by the suffix (0), they 
are connected with the corresponding values of the co-ordinates, <£ Q , etc., by 
the equations 

dY 


d?o = 


#0 


©0 = ^ , etc. 
duo 


( 2 ) 


This strained condition is now to be regarded as initial, and displacements 
from it are denoted by ascribing to the co-ordinates the slightly altered 
values (j )o + $&, 0 0 + 86, etc. For the potential energy of strain we have 

d 2 Y 


V -Y 0 = S</> + ~ ^ my + Sep. BO '.+..., 

dcf)o d0 o dej) o 2 d<f> o dd 0 


(3) 


which is of the first order of the small quantities 8<f>, etc. But Y — Y 0 is 
not now the whole potential energy. In addition to the potential energy 
of strain we have to include that of the steadily imposed forces, represented 
by the terms 

-4>o80-©o80-.... (4) 


The whole potential energy is thus 

(5) 

regard being paid to (2). The total potential energy, as given by (5), is now 
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1906.] On the Dilatational Stability of the Earth 489 

of the second order in S<£, etc., as is obviously required- by the circumstance 
that the strained condition <£ 0 , etc., is one of equilibrium under the proposed 
forces. The coefficients of stability are cPV/dfa 2 , d 2 V / dfaddo, etc., and they 
may differ finitely from the values which obtained previously to the applica¬ 
tion of the forces <P 0 , etc. 

As an example having an immediate bearing upon the matter in hand, let 
us consider the case of a uniform non-gravitating body originally in a state 
of ease. If a small hydrostatic pressure <3> act upon it, the volume changes 
proportionally, and the ratio gives the “ compressibility ” of the body in this 
condition. Under the action of a finite pressure <I>o the volume may be 
greatly altered, especially if the body be gaseous, but the new condition is 
still one of equilibrium and may be regarded as initial. The compressibility 
now may be quite different from before, but it may be treated in the same 
way as depending upon the small change of volume accompanying the 
imposition of a small additional pressure S4>. 

To those who, while accepting the usual elastic theory for bodies in a state 
of ease, repudiate the application to bodies subject to great hydrostatic 
pressure, I would suggest that liquids and solids, as we know them, are not 
really free from stress. In virtue of cohesional forces, there is every reason 
to believe, the interior of a drop of water is under pressure not insignificant 
even in comparison with those prevailing inside the earth, and the same may 
be said of a piece of steel. 

The conclusion that I draw is that the usual equations may be applied to 
matter in a state of stress, provided that we allow for altered values of the 
elasticities. In general, these elasticities will not only vary from point to 
point, but be asolotropic in character. If, however, we suppose that the body 
is naturally isotropic, and that the imposed stress is everywhere merely a 
hydrostatic pressure, so that by pure expansion a state of ease could be 
attained, the case is much simpler and probably suffices for an approximate 
view of the condition of the earth. But although the initial state is one free 
of shear, we are not to conclude that the rigidity is the same as it would be 
without the imposed pressure. On the contrary, there is much reason to 
think that the rigidity would be increased. If there is any analogy to be 
found in a pile of mutually repellant hard spheres, it will follow that an 
infinite pressure will entail infinite rigidity as well as infinite incompressibility. 

In the original draft of this paper I had supposed that it would be possible 
upon these lines to find another and a more practical basis for Professor 
Jeans' analysis. A correspondence with Professor Love* has, however, 
convinced me that this hope is destined to disappointment, and the remainder 
* To whom l am indebted also for other corrections. 
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of the paper loses accordingly much of the interest which at first I felt for it. 
In Professor Jeans’ theory, if A be the dilatation, so that the altered density 
is p (1—A), U the radial outward displacement, E the potential of a volume- 
distribution of density p A, and a surface-distribution of density — pU, the 
displacements £, r\, f are subject to 


<Pj[ 

^ clt 2 




clA 

dx 


+ y“ 


dE 
dx 5 


( 6 ) 


and two similar equations relating to rj, £ In (6) A. and pu are the elastic 
constants of Lame’s notation, and they relate to displacements from the 
compressed initial condition. From equations (6) we obtain, as usual, 

P^ = ^+- 2 P)V 2 &-pV 2 V; ( 7 ) 

and by Poisson’s equation 

y 2 E — — 47T7pA, (8) 

7 being the constant of gravitation. Thus 

P = ( x + V 2A + 4tt7P 2 A, (9) 


which is Professor Jeans’ equation.* 

The solution of these equations is developed by Professor Jeans with the 
view of determining at what point instability sets in. Attention is concentrated 
mainly upon the solution of (9) expressed by a spherical function of order 
one, as being that which bears upon the question of the evolution of the 
moon. 

I had intended merely to indicate a somewhat simpler treatment, following 
more closely the notation and method of Lamb’s memoir, “ On the Vibrations 
of an Elastic Sphere; ”f but as the results so obtained do not agree with those 
of Jeans, it appears necessary to set forth the argument in fuller detail, so as 
to facilitate criticism. 

If in (9) we assume that A is proportional to cos pt, we get 

(V 2 + /^ 2 ) A = 0, (10) 

where w = pM±^ 1- (H) 

and the solution of (10), subject to the condition of finiteness at the 
centre, is 

A = (hr)~ h J n+i (hr ). S n . cos pt, (12) 

J being the symbol of Bessel’s functions, and S n a spherical surface function 
* Loc. cit., p. 162. 

t 4 Proceedings London Mathematical Society,’ vol. 13, p. 192, 1882. 
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of order n. As is well known, J n+h is expressible in finite terms ; in the case 
of n = 0 , (hr) h J k (hr) may be replaced in ( 12 ) by sin hr/hr, a constant factor 
being disregarded. 

Before going further it may be well to consider the particular case of a 
fluid for which yu = 0. Here the solution for A already given suffices to solve 
the problem, and the condition of no pressure at the surface (r = a) gives 
at once 


J n+$(fla) — o, 


(13) 


which with (11) determines p 2 in terms 7 , p, a and the elastic constant X. The 
criterion of stability follows by setting p = 0. In the case of n = 0, where 
the displacements are symmetrical, ha = mir, m being an integer; and we see 
that equilibrium is unstable for symmetrical displacements if 


a 2 p 2 7 > rX. 


(14) 


In general, by ( 8 ) and (10) 


v 2 E = 4 -||£ v 2 a, 


(15) 


so that 



(16) 


where e satisfies throughout the sphere 

V 2 e = 0. 

Substituting the value of E in ( 6 ), we get with regard to (11) 


(17) 



(18) 


k 2 = p 2 p / fjb. 


where 


(19) 


Equation (18) and its companions may be treated as in Lamb’s classical 
paper. A solution is 


1 dA . 1 de 


( 20 ) 


h 2 dx p 2 dx ’ 


where A satisfies (10). In virtue of (17) these values satisfy the relation 




dx dy dz 


and the solution may be completed by the addition of terms u, v, w, satisfying 
(V 2 + ^ 2 ) u = 9, etc., as well as the relation 


du dv dw 
dx dy dz 


Professor Lamb gives the general values of u, v, w. For our present 
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purpose, and with limitation to one order of spherical harmonics, it suffices 
to take 


u = tyn-! (kr) 




J c 2 r 2n+S c l 0 

- 7 ir n+ 1 (kr) - 


( 21 ) 


dx + l 2?i + l .2%+ 3 T "' riV ' dx r 2w+1 9 

and two similar equations, where <j) n is a solid harmonic of degree n ; 
r = y (^ 2 + ?/ 2 + ^ 2 ); and is defined by the equation 


tyn ( 6 ) = 1 * 


2.2w + 3 + 2.4.2?i + 3.2m+5 


( 22 ) 


Save as to a constant multiplier 6 n yjr n (d) is identical with 0“*J w+i (0), as 
employed in (12). \jr n is thus associated with solid in place of surface 
harmonics. The function possesses the following properties 


fa (0) - 2 n+3^ n+1 ^ ; 

(23) 

fa (0) + 2n+l^ n ^ = ^ n ~ 1 ^ ’ 

(24) 

fa(0) fa-i(0) = 2n + 1 '2n + 3+* +l(0) - 

(25) 

A formula in spherical harmonics frequently required is 


r ± _ r2 fd’4> n 2 n+\d 4>n 

9n ~~ 2n + l \dx dxr 2n+1 J * 

(26) 

The term of the nth. order in A is thus 


A w — y]r n (Jir) . (o n 

(27) 

and corresponding thereto 

£_ _ 1 dA n f, de n 

h 2 dx p 2 dx 

(28) 

where u is defined as above, and e m as well as cj) n and co n , is a 

solid harmonic 


of degree n. 

The formation of the boundary conditions to be satisfied at the free surface 
of the sphere (r = a) proceeds almost exactly as in Lamb’s investigation 
(p. 199), the only difference arising from the fact that h has now a different 
value. The first of the three symmetrical surface conditions may be written 




d 


\ dr 


d 


-xk n + [r~— + + = 0. 


(29 


The terms in (29) depending on the parts of £, rj, f which involve A w are 
found to be 

A day | T> d co n 
n dx n df r 2n+1 ’ 


( 30 ) 















Downloaded from rspa.royalsocietypublishing.org 


1906.] 

where 


On the Dilatational Stability of the Earth. 


a _ *■ +V a^nQia) 2n—2. „ ,. 

An ~ a 2^+r h? Vn-lW, 


B n = 


fi 2 n + 

_ _A + 2fi a 2 ”* 3 


, /7 N 2 (»+ 2 ) a 2n+3 . ,, \ 

2 n + 1** (Aa) + 2«+1.2»+3 ^” +1 (/) ' 


In like manner Lamb finds for the terms in (29) arising from u, v, w, 


p d(f) n i-n d (j) n 

n dx + n dx r 2 ” +1 ’ 


where 


C B = -i^Ti ^n(ka)-2(n-l )5 


3\ = 


^ 2 a 2w+3 r 


(*a) + } 'Tea fkd ) j-. 


493 

(31) 

(32) 

(33) 

(34) 

(35) 


n +1 2 % + l 

We have now further an additional part arising from e n , which, it should he 
observed, makes no contribution to A n . In this 

d -I \ de n _ / de n 


dr J dx 


de n . de n d&n 
x ~r~ 4- y - 7 - +z--Z = ne n ; 
dx dy dz 


so that the additional part is 


2n—2 de n 
jo 2 dx 


(36) 


The two most important cases where n = 0 and n — 1 are also especially 
simple, in that (36) disappears. It will be convenient to consider them first. 
When n = 0 ,u, v, w vanish: also, since e 0 is constant, (28) reduces to 


* = 


_ 1 _ dAo 

h 2 dx ’ 


v = — 


1 dA 0 


? = 


1 dAo 
~h 2 dz ’ 


(37) 


h 2 dy 

where A 0 is proportional to yfr 0 (hr). The motion is everywhere purely 
radial. Exactly as in Lamb’s investigation of vibrations without gravity, the 
expression (30) reduces to 

-p d coo 

B °dxT* 

where g> 0 is a constant, so that the surface conditions yield simply B 0 = 0, or 
from (32) 

^±2^ (ha) = ffa (ha). (38) 

Writing 0 for ha, and for yfro and ^ (= — S0~ 1 ^jro / ) their values, we get 

40 


tan 0 = 


4~0 2 


\ -f- 2 [JL 


(39) 
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<r = 0. 

a = -4. 

« = h- 

0-6626 

0-8160 

0-8500 


Except for a slight difference of notation, this is the same as Lamb's 
equation, and his results are therefore available. They are expressed by 
means of Poisson's elastic constant a, and they exhibit Im/tt as dependent 
on a and on the order of the root. To adapt them it is only necessary to 
remember that A 2 , as given by (11), has here a different value from that which 
obtains when there is no gravitation ('7 == 0). On the other hand, although 
7 be finite, hajir may still be equated to Ti/r, where T x is the time occupied 
by a plane wave of longitudinal vibration in traversing a space equal to the 
diameter of the sphere, and r denotes the time of complete oscillation. The 
following are the smallest values of 7 ^/tt corresponding to selected values of 
a , as given by Lamb :— 

O' ~ 

0-8733 

Eor example, if a = \ (Poisson's value), the criterion of stability is 

47T7p 2 a 2 < / 0 . 8160 \ 2 
\ -j- 2 fjb 

If or — 1 the material is incompressible, and motion of the kind now under 
contemplation is excluded. 

When n = 1, (36) again vanishes, though for a different reason from before.* 
The form of the solution is accordingly the same as if there were no 
gravitation. We have from (31), (32) 

* a 2 yjri(ha) 

--' -o- 5 

fjb 6 

B — 2/jl a 5 yjri (ha) 2 a 5 ^ (ha) . 
fJb 3 0 

and from (34), (35) omitting some common factors which have no effect, 

Ci = f 1 (Tea); (42) 


) + p^5 ka ^ 1 ' (* a ) 


The surface conditions (29) are of the form 
dco 1 , X, d co 1 , n d<p 1 


(40) 

(41) 


!ph(fe)-f (43) 


+ “ 1 + C 1 ^ 1 + D 1 A^=o, 

dx ax r 6 dx dx r 6 


and, as Professor Lamb shows, they require that 

Ami + Ci^i = 0; Bi&u -pDi^i = 0. (44) 

It follows that <pi and 001 must be of the same form, and also that 

Ex/Ax = D 1 /C 1 , (45) 

* The terms dejdx, etc., in (28) denote in this case a uniform displacement, as of a rigid 
body, and naturally contribute nothing to the surface condition. 
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in which the values of Ai, Bi, Ci, I)i are to he substituted from (40), (41), 
(42), (43). We find 


^2 (ha) _ 5 2 (^) . 

X + 2 fi (ha) tyi (ha )’ 

or if, in accordance with (23), we replace \jr 2 (6) by —5 6~ l ^i (6\ 


(46) 


_ 4-fA yfri (ha) 
X-f - 2ifj, ha . 'xfri (ha) 


„ x i ■tki. (Ml 

ha . (ha) * 


(47) 


Except for the different value of h , this agrees with Professor Lamb’s 
equation (87). 

In equation (46) there is no limitation upon the value of p. If to find the 
criterion of stability we put p or h equal to zero, 


yjr 2 (ha) = fa (ha) = 1, 

and the equation reduces to 

f 2 (ha) = ^ (Aa). (48) 

4/i 

The equation may also be written in terms of the Bessel functions. The 
relation between J and is 

Ji (x) x fa\(^7 tx) = x fa (x) = sin x, 

J| (x) x (\ttx) — ^x 2 fa(x) = ir" -1 sin—■ cos x. 

Ja (x) x \Z(^ttx) = tV^' 3 ^ (#) = (3(r~ 2 — 1 ) sin aj—3 a" 1 cos £; 
so that in terms of J’s (48) becomes 

(*) = g J, (,0 ; (49) 

or, if we introduce the circular functions, 


(^ 2 ) tana?—3a? 
(T 2 tan (r— 


3_ 

2 o' 


X + 2[jl 

n 


(50) 


Unfortunately (49) does not agree with the result given by Professor 
Jeans. In his notation, when n = 1, 

?/i (from 59) = — -ff , y 2 (from 60) = + T fo '■ 

^ fi 

and (54), (57), (58) give 

<») = W6 J, (*). (51) 

n 2sx 2 


The comparison of processes is rendered difficult by the occurrence of 
several errors (possibly misprints) in Professor Jeans’ paper. Thus (33) does 
not seem correct, and (41), (42) do not follow from (38), (39). Starting 
from Professor Jeans’ equations just mentioned and making use of his (30), 
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(43), (44), (45), (48), (49), I have obtained a result in harmony with my 

equation (49). 

From (50) it is easy to calculate the value of \fp corresponding to any 

value of x. When p — 0, tan x = x, of which the first root is 

0 = 1-4303tt = 4*493. 


This gives an angle of 77 (4-180°). Calculating for angles of 60°, 50°, 

40°, we find 


\f p 

x .. 


00 

3-840 

2-056 

4-493 

4-189 

4-014 


1*221 

3*840 


It seems that the value of x is not very sensitive to variation of p, and for 
such values of the ratio of \/p as are likely to occur, especially under high 
pressure, we might almost content ourselves with the fluid solution (p — 0). 

The simplicity of the cases so far considered, viz., n = 0 and n = 1, 
depends upon our having escaped the necessity of determining the value 
of e n . For values of n greater than unity this function remains in the 
equations, which now demand a more elaborate treatment. From (20), 
( 21 ) 


iL | 1_ d6 n . 

h 2 dx p 2 dx 


/52) 


in which the second term becomes infinite when p = 0. In order to balance 
this, cf> n in (21) must be made infinite of the order p~ 2 or Jc~ 2 . Thus writing 
h~Kt> n = <f > n , we have 

1 dc n _ p dc n i 1_ I" i_ Jc 2 r 2 d^ n 

p 2 dx pk 2 dx k 2 \ 2. 2n 4-1J dx 

__ n r 2n+3 ^ ;2 ' /i2 I \ d 

n-\-l 2714-1.2^4-3 L 2.2?i4-5 ’* J dxr 2n+l ' 


Thus, as in the theory of differential equations with equal roots, we have 
when h 2 = 0, 

1 _ d JL« + u — r2 ^ 4 .P 71 _ T 2n+ * __ d_ e n 

p 2 dx dx p 2.2^4-1 dx p ?^4-l 2n 4-1.2^4-3 dx r 2n+l ’ 


with two similar equations, A denoting again a solid harmonic of degree n. 
From (52), (53) we find for the radial displacement 

u ^Z+yv+zZ _ 1 


dA n nf n p n t & n 
h 2 dr r p 2.2 n 4- 3 


H 

at the surface, where r = a. 


£+■<*“>+ 3 +•"<'“) 


■ r ^fn i P n a s n 

a p 2. 2n 4- 3 


(54) 
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The boundary condition (29) requires a parallel treatment. The terms 
depending on A„ remain as in (30), (31), (32). From (34), (35), we get as 
appropriate for the present purpose 


c ; 

D„ = 


< 2 ”- 2 >{ 1 -2^r x}' 


k 2 a 2 


2 ^ + 1 


-2n — 2 — 


n k 2 o? 
2n+i* 


nk 2 a 2n+s 


~f" 1. 2n -f - 1, 27i ~f" 3 
Equations (33), (36) now give, <I> being written as before for k 2 ^, 


(2n—2) 


pf den_ 1_ d<&, 
fjbk 2 dx k 2 dx 


} ~ no? 

2 ^ + 1 


d& n 


+ 


n a 


,2n+S 


d 


which, when k 2 is made to vanish, is to be replaced by 


dx n + 1.2n + l . 2n+'& dx r 2n+1 ’ 

(55) 


(2n—2)^ + £ — ^ — P 

' /V/V> . . Oai I I rJs\si 


ncv 


2» + 3 


d 


'dx ' fi2n + l dx //, n+1.2n+l. 2n + 2> dx r 2n+1 * 

This is additional to (30). The equations to be satisfied at the surface 


are thus 


A n a, n + (2n-2)f n +P ^~e n . 

fjh 2^+1 




na- 


a»+3 


0 . 


0 . 


(57) 


(58) 


/i n+1.2n+1 . 2n+3 

When n = 1, /„ disappears, and the final condition is found by eliminating 
the ratio a„ : e„ from (57), (58). This would conduct us again to the results 
already arrived at for that case. In general we require another equation 
connecting e„ with m n and /„. 

For this purpose we must recur to the definition (16) of c n and of E. 
A calculation is made by Jeans on the basis of the expression (12) of A 
by means of Bessel’s functions. We have at the surface 

{-U-^ + i)A-i(Aa)-IJ B+i (M).S„-7^i(Aa)-*J/ +5 (Aa).S B ).. 

(59)* 

In order to express this in our present notation (by means of yfr’ s), we see 
by comparison of (12) and (27) that 

(lict) . S w = ( hcdf yjr n (lid) . o) n) ( 60 ) 

so that with use of (24) 

- = (M) 


* In Professor Jeans’ equations (5), (23) the sign of U is positive, but this appears to 
be an error. 
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Eliminating U between this and (54), we find 


fn 


pa% 


(2n+l)p _ 1 _ \ 

H L n(^k + 2pb)h 2 a 2 2.2w + 3J 

2 71 -J-1 




+ 


h 2 


{■ 


tyn-l + tyn 


h 2 ci 2 


, 2n + 3 


'ty'n 




(62) 


The substitution for /„, in (57) now gives 

r* , 2^—2 f 2n + l , , , h 2 a 2 , ~] ' 

n l K + ~W~\ ~ + M .2^ + 3^ +1 J . 

(2n—2) (2ft,+ 1) fi 2n—2 


+ 


pa 2 e n 


ft, (\+ 2 pu) h 2 a 2 2.2 n + 3 ' 2 n +1J 

This equation and (58) determine two values of co n : e n , and the elimination 
of this ratio gives the required final result. We will write (63) for brevity as 


¥a 2 co n + G (pa 2 jp) e n = 0, 
where by (31), and reduction with the aid of (25), 


"F — I 2^—2 

~~ pu 2 ft, +1 h 2 a 2 


{- 


2ft,+ 1 

n 




h 2 a 2 


c _ (2tt-2)(2tt+l)j* _ 2ft—2 

— n(X-\-2fi) h 2 a 2 2.2^ + 3 2?^+l 

Similarly, if (58) be written 

H<x 2 <% + K (pa 2 /pu)e n = 0 , 


2ft,+1.2^ + 3 

n 




we may take 


H = 


X+ 2 p, tyn 


+ 5 


K = 

and the final result is 


2 + 2 ) 

pu 2ft, +1 ' 2ft, +1.2ft, + 3 

n ■ . 


'ifr. 


n+l 


n+1. 2ft,+ 1.2ft, + 3’ 

FK-GH = 0, 


(64) 

(65) 

( 66 ) 

(67) 

( 68 ) 

(69) 

(70) 


giving the ratio yfr n +i (ha) : ^ n (ha) in terms of n, ha, and (X + 2p)!pb. 

In applying results of calculation based upon the assumption of a uniform 
compressibility to the case of the earth where the variation is likely to be 
very considerable, we must have regard to the character of the function ( 12 ) 
by which the dilatation is expressed. When n = 1 or a greater number, (12) 
vanishes at the centre and (when pb = 0) at the surface. The values to be 
ascribed to the elasticities are those proper to an intermediate position, such 
as half-way between the centre and the surface. For a more complete 
treatment we might calculate the balance of the elastic and gravitational 
potential energies on the basis of a displacement still following the same law 
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as has been found to apply to a uniform sphere. In accordance with a general 
principle the result, so calculated, will be correct as far as the first powers of 
the variations from uniformity. 

Another question, interesting to geologists, upon which our results have a 
bearing is as to the effect of denudation in altering the surface level. The 
immediate effect of the removal of material is, of course, to lower the level, 
but if the material removed is heavy and the substratum very compressible, 
the springing up of the foundation may more than neutralize the first effect 
and leave the new surface higher than the old one. So far as I am aware 
discussions have been based upon the elastic quality merely of the interior 
without regard to self-gravitation; but, as is easy to see, if the condition be 
one approaching instability, the effect of a pressure applied to the surface may 
be immensely increased. 
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